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■ We study antiferromagnetic spin chains with unfrustrated long-range interactions that decays as 

a power law with exponent /3, using the spin wave approximation. We find for sufficiently large 
spin S, the Neel order is stable at T = for j3 < 3, and survive up to a finite Neel temperature for 
(5 < 2, validating the spin- wave approach in these regimes. We estimate the critical values of S and 
T for the Neel order to be stable. The spin wave spectra are found to be gapless but have non-linear 
momentum dependence at long wave length, which is responsible for the suppression of quantum 
and thermal fluctuations and stabilizing the Neel state. We also show that for /3 < 1 and for a large 
but finite-size system size L, the excitation gap of the system approaches zero slower than I/ -1 , a 
behavior that is in contrast to the Lieb-Schulz-Mattis theorem. 
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I. INTRODUCTION 



Antiferromagnetic (AF) spin chains have attracted considerable interest from physicists in the last two decades, and 
continue to be subjects of active research at present. There are several reasons why they are of such strong interest. 
Firstly, quantum antiferromagnetic spin chains are important examples of a larger class of strongly correlated systems, 
whose ground state and low-energy behavior differ from their higher dimensional counterparts in qualitative ways. In 
the case of AF spin chains, quantum fluctuations destroy the Neel order in the ground state no matter how big the 
I \ the size of the spin is, while in higher dimensions the Neel order is stable regardless of spin size, in the absence of 
frustration. Secondly, the spin chains are of interest to physicists because they are ideal playgrounds for various types 
of theoretical approaches. A prominent example here is the work of Haldane 1 , who mapped the AF spin chains to 
quantum non-linear sigma models, and predicted that the integer chains have a gap in their excitation spectra while 
no gap exists for half integer chains, based on the absence or presence of a topological term in the mapping. This 
fundamental difference is consistent with, and to certain degree implied in the Lieb-Schultz-Mattis (LSM) theorem-*, 
which states that for a Heisenberg AF chains with length L and periodic condition, for half-integer spins, there exists 
an excited state with energy separated from the ground state that is of order l/L; no such theorem exists for integer 
chains however. 
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' The studies of AF spin chains, and the results mentioned above, are restricted to models with short-range in- 
(SJ ■ teractions. In this work we study AF chains with interactions that decay as power laws and without frustration: 

o 

p : h = X : | !< ■'■<..,». ■ sj, (i) 



with 



JiJ = J/\i~jf, (2) 

where J > determines the overall energy scale of the system and (3 is the power-law exponent that controls the 
O ■ decay of the interaction. The factor (— ensures that spins sitting on opposite sublattices have antiferromagnetic 
interactions and those sitting on the same sublattice have ferromagnetic interactions, thus there is no frustration. Our 
motivation of the study comes from the following considerations. Firstly, such power-law long-range interactions can 
in principle be realized in experimental systems; one example of which being the RKKY 3 interaction mediated by 
conduction electrons that decay as power laws, with an exponent that depends on the details of the conduction 
electron Fermi surface. Secondly, as we will show, such long-range interactions tend to suppress quantum as well as 
thermal fluctuations, thus increasing the range of interaction has an effect that is somewhat similar to increasing the 
dimensionality of the system. On the other hand the dimensionality is discrete while the power-law exponent for the 
interaction can be tuned continuously, thus providing a tuning parameter for the fluctuations; it is of interest to study 
how the system behave under such tuning. 

Anticipating the stability of Neel order in the presence of such long-range interactions, we study the models using 
the spin-wave method. We obtain the following results, (i) We show the Neel order is stable at zero temperature for 
(3 < 3 and sufficiently large S, justifying the usage of spin- wave method in this case. We also estimate the critical size 
of the spin for the Neel order to be stable, as a function of (3. (ii) In this case the spin-wave excitation spectra take 
the form Uk ~ k 1 in the long wave-length, with 7 < 1 and varying continuously with (3. (iii) Extending the spin-wave 
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calculation to finite temperature, we show that the Neel transition temperature Tjy is zero for (3 > 2 while finite for 
(3 < 2. We determine Tjv as a function of S and /3. (iv) For a finite-size system with size L and periodic boundary 
condition, and (3 < 1, we find the lowest excitation energy approaches zero slower than 1/L as L increases for both 
half-integer and integer spins, thus "violating" the LSM theorem. Of course the LSM theorem applies to spin chains 
with short-range interaction only; here we have provided explicit examples of how it is invalidated by the presence of 
long-range interaction. 

The remainder of the paper is organized as follows. In Sec. II we discuss the application of spin wave technique 
to this model. In Sec. Ill and IV we present and discuss the significance of our results. In Sec. V we summarize our 
work and discuss the implications of our results. 

II. THE SPIN WAVE APPROACH 

We consider a Heisenberg antiferromagnetic chain with unfrustrated power-law long-range interaction with the 
Hamiltonian given by Eq. The central issue we address in this work is the stability of Neel state at zero or 

low temperature. It is thus natural to use the spin-wave method based on the Holstein-Primakoff transformation^ 
that maps spin operators to boson operators, and check its self-consistency. The procedure is rather standard^ we 
nevertheless include the details here for the sake of completeness and establish notation for later treatment. We 
divide the chain into two sublattices and represent the spin operators in terms two types of bosons : a bosons which 
live on A sublattice and and b bosons which live on B sublattice. Up to order 1/5, where S is the size of spin, the 
Holstein-Primakoff transformation for the spin operators can be written as the following : 

S? =S- a f a, Sf = \/2Sat(l - a t a/(2S')) 1/2 ~ V2Sol; i € odd 

S? =-S+ bh, = V25(l - tfb/{2S)) 1/2 b ~ V2Sb; i E even. (3) 
Using this transformation, the Hamiltonian in Eq. can be separated into three terms as follows : 

H Hodd—even H-odd—odd H-even—evem 

where H dd-even, H oc id-odd, and B dd-even are defined as : 

L 

Hodd-even = J]] ^2i-l,2j [~ S 2 + S (a\ i _ l a2i-l + b\^>2j + 0,2i~lb2j + a 2i-1^2j)) 

i,j 

L 

Hodd-odd — — J]] J2i-l,2j-l\—S 2 + S{a\ i _ 1 a2i-\ + a 2 j_ 1 d2j-l + &2i-\a\j-\ + a 2i- 1 a 2j-l )] 
i<j 
L 

H eve n-even = ~ ^ hi,2i \S 2 + S(b' 2i b 2 i + 4^2j + b 2 ib\ j + 4i & 2j)]' ( 4 ) 

i<j 

We diagonalize this quadratic Hamiltonian by going to momentum space and then diagonalizing by a Bogoliubov 
transformation: 

H = constant + JS ^ [(a - f{k))(a\a k + b\b k ) + g(k){a\b\ + a k b k ) (5) 



where 



L/2 

a = 2 lim > r-r, 

n=l y ' 
L/2 

f(k) = 4 lim ^ [cos(2nfc) - 1], 

n— 1 v ' 

g (k) = 2^g^l IF cos(2n-l) fe ; (6) 
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using the Bogoliubov transformation, the Hamiltonian (JSJ can be diagonalized and be written in terms of free boson 
operators Cf. and d^: 



H = constant + JS'S^uik(cj l Ck + d' k dk 



(7) 



where 



= V(a-/(fc))2-( 5 (A)) 2 . 
The correction to staggered magnetization is given by : 

Am = — < a l a k >— Am, + Am T (T), 



(8) 



(9) 



where Am, and Amr (T) , which represent the quantum and thermal fluctuation correction respectively, are given by 



Am„ 



dk 1 
2^2 



a-f(k) 



. dk 

Am T (T) = / - 



CO). 

a - f(k) 



1 



oE k /k B T _ l ■ 



(10) 



We will visit these equations frequently when we discuss the validity of the spin wave approach later in the text. 

It is clear that the correction to magnetization is dominated by the small k behavior of the spin wave spectrum. 
We thus need to obtain the small k behavior of the expressions given in Eq. © . To do that we express them in terms 
of the Bosc-Einstein integral function^ defined as : 



F(a,v) 



1 



„a-l 



dx- 



r(a) J ™ e x+v - 1 l a ' 2 a 
and rewrite the cos(nfc) term in f(k) and g{k) as the following : 

cos(nfc) 



E 



E 



(ii) 



(12) 



The analytical properties of F(a, v) near v — are known and are given by : 
F(a,v) = r(l - a^"- 1 + V C( " " n) (-v) n , (a d Z) 



F(a,v) 



(-«) 



Tyr[Ej-^H + E 



(a-l) 



C(a - n) 



(-v) n ,(aeZ) 



(13) 



n^a— 1 

where C(s) is the zeta function. We will use these properties in our later treatment. 



III. SPIN WAVE SPECTRA AND CORRECTIONS TO STAGGERED MAGNETIZATION 

In this section we analyze Eq. @ for different values of j3 to obtain the spin wave spectra and calculate the correction 
to staggered magnetization, to determine the validity of the spin wave approach. 



A. p > 3 

Equations 1|12|) and (|13f) are the main ingredients to analyze Equation © which can be summed up in closed forms. 
Up to leading order in k the relations in Eq.@) for (3 > 3 read: 



* (2n - \y 



2(l-2-")CG8) 



4 



/(*) = 4 



\ " 



1 



71=1 V ' 

2 2 "^ [R(-F(/3 
cfc 2 



cos(2?ifc) — 1] 
2**)) - CG9) 



^ (2n- l)/ 3 

71=1 V ' 

rcos(nfc) 



cos(2?i — l)fc 
cos(2rtfc) 



(2n)f 



a 



c'k 2 , 



(14) 



where c and d are positive constants. The same results can also be obtained by expanding the cos(nfc) term to order 
fc 2 in/(fc): 



y~\ cos(nfc) 



1 



n 







(15) 



it 



in which the sum converges as long as /3 > 3; together with a similar expansion for one reproduces Eq. I|14[l. 
The spin wave spectrum can be easily shown to be linear in k: ujk oc fc, and the T = correction to the staggered 
magnetization from long-wave length spin wave fluctuation: 

f dk , 
Am, ~ / — (16) 
J u k 

diverges logarithmically for f3 > 3. This immediately indicates that the spin wave approach is not valid for > 3 at 
zero temperature. The results obtained here are essentially the same as the spin wave calculation for nearest-neighbor 
interactions only£ 

For j3 = 3 the expansion we did above is no longer valid because the sum is divergent. We rely instead on the 
Bose-Einstein integral function as defined in Eq. I|llfl to calculate Uk and Am q . After a little algebra we find 
(jJk ~ fey] log(fc) | which leads to the correction of staggered magnetization that diverges as y/\ log(L)|, where L is the 
system size. We thus conclude that the quantum fluctuation destroys Neel order, and the spin wave approach is not 
valid for (3 > 3. 



B. 



1< P < 3 



We now turn our attention to the case 1 < j3 < 3. As in (3 = 3 case we are no longer able to expand the cos(nfc) 
term in f(k) and g(k) because the coefficient of k 2 is divergent so we again take advantage on the mapping onto the 
Bose-Einstein integral function. In the long wave-length regime, the relations given in Eq. (JSJ read : 



;(2n-l)* 



2(1 - 2^)C(/3) 



/(*) 



9{k) 



OO 

71=1 V ' 

2 2 -^{F((3 1 -2 l k))~m 

i 



OO 



cos(2n — l)fe 



^ (2n- lY 

n=l y > 

rcos(nfc) cos(2nfc) 



(2n)/3 



(17) 
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where the function (f>(f3) is given by : 



i 



T(f3) cos[7r(/3-2)/2]' 

with r(/3) being the gamma function. The long wave-length spin wave spectrum is given by : 



which is sublinear, and T = correction to staggered magnetization by : 



Am„ ~ — — 



1 

2n 



a(f3) 



3</>(/3) 3-/3 



r (3-/3)/2 



3a(/3) 1 + /3 



7r (l+/3)/2 _ ^ 



(18) 



(19) 



(20) 



which is convergent for /3 < 3. These results show that the system supports gapless excitations, the spectrum follows 
a sublinear power-law at small momentum k, and that the Neel order at zero temperature survives, for large enough 
S, for 1 < j3 < 3. Our results agree with an earlier work presented by Parreira, Bolina, and Pereaiwho show the 
existence of Neel order for f3 > 3 and the presence of Neel order for < 3 at zero temperature using rigorous proof. 
However the excitation spectra were not studied in this work, nor was the critical value of 5* for the stability for Neel 
order calculated. Another support for our results at zero temperature is offered by the work of Aoki& who studied 
the same model we are studying for the case (3 = 2 in ID and 2D using spin wave theory In that work he found that 
there exists Neel order at zero temperature in 1 dimension for f3 = 2 which is in agreement with our conclusion. 

We may also estimate the critical size of the spin, S c , above which the Neel order survives, by setting the correction 
to the staggered magnetization equal to the spin size: Am, = S c . As (3 — > 3, Am, is dominated by long-wave length 
spin- wave fluctuations, and we obtain 



3(f)(0) 3-/3 



^(3-/3)72 + 



3a(l3) 1 + (3 



7r d+0)/2_ 7r 



0.41 

V3 - T 



(21) 



a result we expect to be asymptotically exact in the limit (3 — ► 3. On the other hand we also find the quantum 
correction gets suppressed very rapidly as /3 decreases from 3; for examples we find S c — 1/2, for (3 = 2.63 and S c — 1 
for (3 = 2.85, suggesting the Neel order would survive for any spin for (3 < 2.6. 

We also calculate the correction to staggered magnetization at finite temperature. First we discuss the case for 
[3 > 2. The thermal correction to staggered magnetization is given by : 



Ato t (T) 



k B T 
2irJS 



dk 



30(/3) 



1 

3a 



(22) 



which diverges as L^~ 2 for (3 > 2. For (3 = 2, it is a simple exercise to show that the spectrum behaves like uju ~ Vk, 
and the small k contribution to the thermal correction of staggered magnetization diverges as y/\ log(L)|. These results 
indicate that thermal fluctuations destroy the Neel order for (3 > 2 at any finite temperature. They are consistent 
with an extension of the Mermin- Wagner theorem that Bruno advanced, 9 which proves the absence of Neel order 
for (3 > 2. For classical antiferromagnets in ID, it has been shown, using Monte Carlo simulation, that there is no 
magnetic ordering at finite temperature^ 

For (3 < 2 the correction to staggered magnetization is given by : 



Am T (T) 



k B T 
irJS 



r(2-/3) 



3(2-/3)0(/3) 



7T 

3a 



(23) 



This convergent correction shows that the Neel order survives at finite temperature for (3 < 2. The Neel transition 
temperature Tjv can also be estimated by applying the same rationale used to estimate the critical value of S at zero 
temperature. By using Eq. I|23|l we find : 



T N (S,/3) 



In the limit (3 — > 2, we find Tjv vanishes linearly: 



irJS 



r 2-/3 



3(2 - I3)(f)((3) 3a 



3ir 2 JS 
ks 



(2-/3) 



(24) 



(25) 
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We see that increasing the range of interactions (or decreasing (3) in the chains has effects that are similar to 
increasing the dimensionality of the systems. For (3 > 3 we find absence of Neel order at both zero and finite 
temperature, a genuine one-dimensional (ID) behavior. For 2 < (3 < 3 we have finite Neel order at zero temperature 
which gets destroyed at any finite temperature, similar to the 2D situation. Finally for (3 < 2 the Neel order is stable 
at zero and low-enough finite temperature, a behavior expected for dimensions above two. 

In contrast to the antiferromagnetic case we are studying here, the ferromagnetic models with long range interactions 
have been studied more extensively. Classical Heisenberg model with long range ferromagnetic interactions has a 
phase transition at finite temperature in I dimension when the interactions decay slower than 1/r 2 . There is no phase 
transition at finite temperature when the interactions decay faster than l/r 2 iii This result for the classical case in 
1 dimension is confirmed by Monte Carlo simulation^ The quantum Heisenberg model with long range interactions 
have also been studied using the modified spin wave theoryii&ii It was shown that there exists a magnetic ordering 
in 1 dimension as long as the interactions decay slower than I /r 2 . 



In this section we consider the case < 1. The reason we separate (3 < 1 case with the rest is because there are 
divergences in the thermodynamic limit which require special care in their analysis. Physically, this is closely related to 
the fact that the ground state energy grows faster than the system size (i.e., it becomes "superextensive" ) , if the local 
energy scale J is not rescaled according to the system size. For this reason we will not discuss the finite temperature 
(or thermodynamic) properties of the system, as the definition of temperature becomes somewhat ambiguous; we will 
focus instead on the ground state properties of the system, which is free of such ambiguity. 

For the reasons mentioned we need to work explicitly with a finite system size L, defined as the number of spins 
per sublattice (so the total number of spins is 2i), and treat k and L as two independent variables. For a start, the 
summation in a 



C. /3 < 1 



L/2 



I 




(26) 



diverges for [3 < 1 if we run the summation to infinity. For large but finite L, we have 




Similarly 




(27) 




2r(I - 13) cos(tt(/? - l)/2)kP~ 1 - 2L 1 -^/{\ -(3) /3 < 1; 



13=1, 



and 




(28) 




The spin wave spectrum reads : 




(29) 



7 



E{k) 



3J51og(Z)(l + 61og(fc)) (3 = 1 
3JS'L 1 -' 3 (1 - fefc^- 1 ) /3<1 



which approach L-dependent constants as fc — * 0. Here 6 oc l/log(L) for (3=1 and b <x L 13 1 for /? < 1. Correction 
to staggered magnetization at zero temperature can be calculated easily using the relations derived above to yield : 



Am. 



1 



" log(L) 
1 



(3=1 
(3<1, 



(30) 



suggesting the quantum fluctuation gets completely suppressed as system size grows. 

For (3 = the calculation becomes particularly simple; the relations for a, f(k), and g(k) in Eq. © become : 



a = L, 



/(*) = E [■ 

9{k) 



e ik S 2 + e -ik S 2 _ 2 



2i(4,o - 1), 



Si 

= L5k,o- 

The spin wave spectrum for k ^ is given by : 

E k = JSLyJ{l + 2) 2 = 3JSL, 
which is k- independent, and the correction to staggered magnetization is given by 

Am„ 



We will compare these with an exact solution for this special case in the next section. 



(31) 



(32) 



(33) 



D. (3 = : exact solution 

The infinite range (/? = 0) antiferromagnetic chain with no frustration is given by the following Hamiltonian : 



2L 



H = J^2(-l)*-H +1 Si ■ S 

ij 

which can be solved exactly in the following manner. We introduce: 

Sa = ^ Si, 

sb — Si, 



(34) 



ieB 



where Sa(Sb) is the total spin operator for sublattice A(B), to rewrite the Hamiltonian in the following form 



H = J 



S A ■ S B - ((S 2 A + S|)) + ( J](^) 2 + ]T(S 4 ) : 



ieA ieB 

We define the total spin operator Stot — S a + Sb to further simplify the Hamiltonian given above to become 



H = J 



S 2 tot -l(S 2 A + S%) + 3L/2 



(35) 



(36) 



(37) 
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The Hamiltonian in Eq. I|37|) can be diagonalized in the total- £> basis of states given by \(Sa, Sb)', Stot > where 
Sa(Sb) and S to t are the total spin quantum number in sublattice A(B) and in the system respectively. Using this 
basis, the energy can be easily obtained as : 



E = J 



\s to t{Stot + 1) - \{Sa(S a + 1) + S B (S B + 1)) + 3L/2 . CM 



To minimize the energy we must have all spins aligned in each sublattice and have a minimum of Stot- This means 
that Stot = and Sa = Sb = LS, where S is the spin size, will minimize the energy and give us the ground state. 
The momentum quantum number of the ground state is (tt) for even (odd) L. The lowest energy excited state is 
obtained by having Stot = 1 while still maintaining maximum Sa and Sb- The energy gap is given by : 

AE = E ex - E gs = J. (39) 

This particular excited state has a momentum quantum number that differs from the ground state by 7r, which 
corresponds to momentum k = in the spin wave approach, due to the doubling of the unit cell in that approach. 
We will say more about this in the next section. To obtain excitations with generic k however, we must change 
either the Sa or Sb quantum numbers. There exist two branches of degenerate low-lying excitations, corresponding 
to Sa — LS — 1 or Sb = LS — 1 and Stot = 1, with excitation energy 

AE = E ex - E gs = J(l + 3LS), (40) 

which grows linearly with system size, and has no fc-dependence. This result agrees with the spin wave solution 
obtained earlier in the limit S — > oo, as expected. 



IV. EXCITATIONS AT k = AND STATUS OF THE LIEB-SCHULTZ-MATTIS THEOREM 

The Lieb-Schultz-Mattis (LSM) theorem^ states that for half-integer spin chains with length L and short-range 
interaction, there exist an excited state whose momentum differs from the ground state by tt, with energy that 
vanishes at least as fast as 1/L as L — > oo, 2 Recently the theorem has been extended to spin chains with power-law 
long range interaction, and it was found that the theorem remains valid for (3 > 2^i& The situation is unclear for 
0< 2. 

In this section we check if the LSM behavior still holds for (3 < 2, using the spin- wave method. As discussed above, 
due to the doubling of the unit cell, excitations whose momenta differ from the ground state by either tt or show 
up as k — excitation in the spin-wave approach. If one blindly use the linear spin-wave results however, one would 
always find E k=0 = 0. But this is an artifact of the linear spin- wave approach which maps the k = modes to 
harmonic oscillators without a restoring force. Thus in order to study the excitation that are relevant to the LSM 
theorem, we must treat the k = modes more carefully. 

To do that, we start by rewriting the Hamiltonian as given in Eq. in the momentum space: 

H = J2 E J (^fc • S B k e-^ ( J &) S £ ■ S- k e- ih S2 + J(S 2 )S* ■ S B k e^ k •*») , (41) 

k Si k 62 

where 

Sf /B = ^fj:si /B e- ik -°», (42) 
V7j k 

and A(B) denotes odd(even) sublattice. Instead of applying the Holstein-Primakoff mapping for all terms in H , we 
separate out the k = term in H and apply Holstein-Primakoff mapping to the k ^ terms only. Since to linear 
order the k = term commutes with the other terms in H, they can be diagonalized independently. The spin wave 
treatment for the k 7^ terms gives the spectra obtained earlier, except that k must be nonzero. On the other hand 
the k = term 

Si i£A i£B S 2 ieA ieA 

= i J(Si)Sa -Sb-jJ2 J ^ ( S a + S b) . ( 43 ) 

<5i S 2 
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takes a form identical to the Hamiltonian for (3 = that was solved exactly in the previous section. We can easily 
solve this Hamiltonian to obtain the excitation energy at momentum 7r measured from the ground state momentum, 
or k = for the doubled unit cell: 

A£=^, (44) 

where a depends on the power law exponent (3 and is given by Eq. For (3 > 1, a is convergent in the large L 
limit and is given by Eq. (|14|l . This means that the energy of the excited state vanishes as 1/L as L — > oo. For (3=1, 
a diverges as ln(L) as shown in Eq. (|26[) and the energy vanishes as \og(L)/L. For j3 < 1, a diverges as L 1 ^ 13 as 
shown again in Eq. I|2tj[l and the excitation energy vanishes as L~@ . We thus find that the LSM behavior holds for 
1 < (3 < 2, despite the the absence of a proof for this range of (3. On the other hand the LSM theorem is "violated" 
for (3 < 1. 



V. SUMMARY 



We have studied antiferromagnetic chain with unfrustrated long range interaction using the spin wave technique. 
We find that this approach is valid for (3 < 3 at zero temperature for sufficiently large size of spin, and (3 < 2 for 
sufficiently low finite temperature, due to the stability of Neel order. Within the range of validity of this approach 
we find that the system has a gapless excitation and the excitation spectrum follows a non trivial k dependence. We 
also study how the excitation gap closes in this system in the limit L — ► oo, and find a behavior that is in contrast to 
that predicted by Lieb-Schultz-Mattis theorem for chains with short range interactions, when [3 < 1 . 
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